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We present a simple and straightforward method for relating the form 
factors in HQET, as defined by the covariant trace formalism, to the overlaps 
of the rest frame wave functions of the light degrees of freedom. We also point 
out several inconsistencies present in recent calculations of the radiative rare 
B decays, and also show how these can be fixed even within the framework 
of the non-relativistic quark model. 



Abstract 



1 Introduction 



It is by now a well established fact that hadronic systems containing a single heavy 
quark (rriQ 3> Aqcd) admit additional symmetries which are not present in the full 
QCD Lagrangian The light degrees of freedom (LDF) in such hadrons typically 
have four-momenta small compared with the heavy quark mass. For these systems it 
is appropriate to adopt an effective theory (HQET) in which the heavy quark mass 
goes to infinity, with its four- velocity fixed [fj, |||] . Since an infinitely massive heavy 
quark does not recoil from the emission and absorption of soft (E w Aqcd) gluons, 
and since magnetic interactions of such a quark are negligible (~ ~)j the strong 
interactions of the heavy quark are independent of its mass and spin. Because of 
this, HQET leads to relations between different form factors describing transitions 
in which a hadron containing a heavy quark Q and moving with four-velocity t> M , 
decays into another hadron containing a heavy quark Q', and moving with four- 
velocity v ,fl . In this way the number of independent form factors for these decays 
is significantly reduced. For example, the six form factors describing semileptonic 
decays B — > D^eVe are, in the heavy quark limit, reduced to a single unknown form 
factor, the Isgur- Wise function (IW) £ (v ■ v') . 

Since the unknown Lorentz invariant form factors describing a particular decay 
cannot be calculated from first principles, one still has to rely on some model of 
strong interactions in order to estimate them. In general these form factors will 
be related to the overlaps of the wave functions of the LDF in the hadrons before 
and after the decay. However, one has to be careful in identifying form factors 
directly with the overlap of the two wave functions. Depending on the definition 
of the particular form factor, there may be a coefficient of proportionality involved. 
If these coefficients are not taken into account significantly incorrect results can be 
obtained no matter which model for the wave functions one uses 

In [|J it was observed that the heavy quark limit implies a simple formula for 
the wave function of any particle containing one heavy quark. Based on this we 
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offer a straightforward method for relating the form factors, as defined within the 
framework of the trace formalism |5], ^ 0], to overlaps of wave functions of the LDF 
before and after the decay. Even though we are interested here only in mesons, it 
is obvious that an analogous calculation can be easily done for baryons. The paper 
is organized as follows: in Section |] we review basic definitions of the covariant 
trace formalism. Section ^ describes how one can easily identify the form factors in 
terms of overlaps of the LDF. We also give results for the several cases of interest 
(transitions of 0~ state into excited states). As a simple application of the results of 
this paper we consider the non-relativistic quark model in Section f|. Among other 
things, we also make a few comments on some recent calculations of rare B decays 
into iT-resonances || || . Our conclusions are summarized in Section 



2 Covariant representation of states 

The counting of the number of independent form factors is most conveniently done 
within the framework of the trace formalism, which was formulated in |5], |5| and 
generalized to excited states in [^]. Following 0, and using notation of [|J, the 
lowest lying mesonic states can be described as follows: 
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G(v,e) = ^^(^ + 1)^-^(7^ + ^)], J P = r,j = l, (7) 

G*(v,e) = ^{jb + ljw"", J p = 2-,j = ^. (8) 

Here m and i> are the mass and the four- velocity of the heavy meson, and j denotes 
the total spin of the LDF. Also, e M is the polarization vector for spin 1 states (sat- 
isfying e • v — 0), while the tensor & v describes a spin 2 object {e^ u = e Ufl , e fiu v u = 
, e M M = 0). These eight states form four doublets: (C, C*) is the L = doublet, 
(£", £"*) and (F, F*) are the two L = 1 doublets, and (G, G*) is the L = 2 doublet. 

Matrix elements of a heavy quark current J(q) = QTQ between the physical 
meson states can be calculated easily by taking the trace (u = v ■ v'), 

(V'(v')\J(q)\V(v)) = Ti[M'(v')TM(v)}M(u) , (9) 

where M' and M denote appropriate matrices from (|I])-(^|), M = ^y°M^ , and 
M.{uS) represents the LDF. Again following 0, ||, we define the IW functions for the 
transitions of a 0~ ground state into an excited state by 

£o(w) , 07 -> (07,17) , 

2 2 2 

£ E (u) , 07-^(01,11) , 



M(u) 



(10) 

07^(lt,2f 



J 3 

2 2 2 



£g(^K , 07 -> (17,27) . 

2 2 2 

The vector index in the last two definitions will be contracted with the one in the 
representations of excited states We now show how one can define the IW 

functions given above in terms of the overlaps of the wave functions of the initial 
and final states of the LDF. 



3 Defining IW functions 

It has been pointed out Q that the assumption of the heavy quark limit implies 
a simple formula for the wave function of any particle containing one very heavy 
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quark (with total angular momentum J and its projection A), 

*5aV) = E U^j^XqIJX^jx^ux^v) . (11) 

In this formula (a) refers to all other quantum numbers of the meson, u\ Q (v) is the 
free Dirac bispinor describing a heavy quark with spin |, helicity Aq, and velocity 
v (and normalized to uu = 2m). is the wave function of the LDF with total 
angular momentum j (in the rest frame of the particle), and its projection Xj. For 
a meson, this is the wave function of the light antiquark. 

From (|TTD it can be easily seen that matrix elements of the heavy quark currents 
J(q) = Q'TQ are linear combinations of matrix elements 

^' jIXj Xv')\^ 3X] {v))u\ qI {v')Tu Xq {v) . (12) 

For a given T, u'Tu is a product of known matrices, and therefore all the unknown 
dynamics is contained in the overlaps of LDF wave functions (^'l*^). We choose 
the spin projection axis (z) as the velocity direction of meson fy' as seen in the rest 
frame of meson From the independence of the overlap on the direction of the 
x-axis we then have 

($;, v I$ jAj ) = 0, if A^A,-. (13) 
Similarly, from the independence of the overlap on the orientation of the y axis, 

W'V l«W = '/'/'( D r a . a,) , (14) 

where r/ and rj are the orbital parities of the initial and final state of the LDF. In 
the case of interest to us, for the change of the orbital angular momentum from L 
to L', rjT]' = (— 1) L+L ' . Equations (|13|) and (|14D have been given in |4j]. 

As an illustrative example, we choose the CT — > + transitions and axial-vector 
current (r = 7^75). From (|TT| ) we have 

^ = ^tU^ ~ ^tKu,} , (15) 

\/2 2 '2 2 2' 2 2 
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where + and — refer to + and states, respectively. Also, from ( pT3[ ) and ([L4]) one 
can see that 

(^ + K(v')\$t\(v)) = -(^\\v')\<S>[~\(v)} , (16) 

2 > 2 2 > 2 2 ' 2 2'2 

and all other overlaps are zero. Therefore, it immediately follows that 

(0V/|r|(T, V > = hu>iTu_i -u>fui}($'tl(v')\$tl(v)) . (17) 

Z 2 2 2 2 2 '2 2'2 



Now choosing T = 7 75 and evaluating (|T7|) in the rest frame of meson, where 
i> M = (1, 0, 0, 0) and v'^ = (lu, 0, 0, y/u> 2 — 1), one easily obtains 

(0 + ,v'\^ 5 \0~,v) = -v / mm'v / 2v / ^TT($ / i + i ) (v , )|$ ( rl (v)} ■ (18) 

2 ' 2 2 ' 2 

On the other hand, using (|j), @ and flip] ) inside @ one finds 



(0 + , w'|7 M 7 5 |0", v) = Vnirn'i-v 11 + v'^ E (u) , (19) 
which specialized to the rest frame of 0~ state yields 

(0 + ,v'\-f s %\0~,v) = Vmm'Vu 2 - l&H . (20) 
Comparing and fl2"0|) we obtain (apart from the irrelevant overall sign) 

= • (21) 

yCJ — 1 2'2 2'2 

Of course, in order to obtain this formula for £e (with the same overall sign), 
we could have chosen any other component and any other current giving a non- 
vanishing matrix element. Also, instead of 0~ — > + transitions, we could have 
chosen 0~ — > 1 + transitions, and specialize to any of the three possible polarizations 
of the 1 + state. Finally, any other reference frame besides the rest frame should 
yield the same expression. 

Let us summarize the results obtained following the simple proceedure outlined 
above for several cases of interest. We emphasize that all the results given here were 
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explicitly verified for the choices of T = 7 M , 7^75, 7^7^, 7^7^75, and in two convenient 
reference frames: besides the rest frame of the 0~ meson, we have also used the Breit 



frame (v = -v'), in which v" = ( ^f 1 , 0, 0, -yj *=±) and 1/" = (^/^, 0, 0, ^f 1 ). 
Polarization vectors describing spin 1 states (e • t/ = 0) were the standard ones, 
e^' = =F^(0, 1, ±i, 0) in both frames, = (\/uj 2 — 1, 0, 0, to) in the rest frame of 



, and = (y 0, 0, y in the Breit frame. Let us also, for the sake of 
simplicity, define 



(^|$> = (CV(t/)|M 0, i(«)> , (22) 
J ' 2 2 > 2 



' 2 

and state our results in terms of this overlap: 

• 07 — > (07, 17) transitions. In this case we obtain 



fcM = y^W&> ■ (23) 
This expression was obtained in [|lOf| . 

07 — > (0t, it) transitions. Here, as shown in the previous section, 



&( W ) = ,/_?_<$'|$) . (24) 



0i — > (it, 27) transitions. In this case we have 



£f{w) = J — ^— -—($'!$) . (25) 



w — 1 a; + 1 



0i — > (13,23) transitions. Here, 



fo(w) = v-rT^ L rW> • ( 26 ) 

V Co> + 1 Cl) — 1 
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We now proceed to relate the overlaps ($'|$) to the actual wave functions de- 
scribing the LDF in the rest frame of the particle. Following |T0| , IT |, in the valence 
quark approximation, we can write for the LDF wave function in the rest frame of 
the particle 

$(°)( x ) = ( /)(°)(x)e- i ^ i , (27) 

where Eg denotes the energy of the LDF. The LDF wave function of a meson moving 
with (ordinary) velocity (3 along the z axis (laboratory frame) is then given by 



$(x') = S{f3)& 0) {x) , 



(28) 



with x' = hr 1 (f3)x being the laboratory frame, x the rest frame of the meson, and 
S(f3) is the wave function Lorentz boost. 

Since IW functions are Lorentz invariant they can be calculated in any frame. 
Particularly convenient is the Breit frame, in which the two mesons move with equal 



and opposite velocities. As noted in [10|, the wave functions relevant for the overlap 
($'|$) are at t' — in the Breit frame. Therefore, denoting the three-velocity of the 



final meson as (3, by the use of (p8|) we have 



($'(t/)l $ 00> = / d 3 x'$ ,j! (x')$(x%, = o 

= J d 3 x'<5>'^(x + )S\(3)S(-(3)<S>(x_)\ t , =0 

= |rfV$ / W t (x + )$(x_)| i , =0 . (29) 

In this expression x + and x_ denote the rest frames of the final (moving in the +z 
direction) and initial meson (moving in the -z direction), respectively. To obtain 
the last equation we have used the fact that Lorentz boosts satisfy S^(f3) = S(f3) = 
S' _1 (— f3), so that boost factors cancel out. Also, we have (/3 = |/3|) 



x±\ t '=o = A(±(3)x' t , =0 = (TlPz^x^y'^z') 



(30) 
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where, in terms of w, 



7 = (31) 



Using (|3) and (0) in (||) we find 

($V)i$oo> = / d 3 x> ,(o)t (x^y / ,7^)0 (o) (^^z/^7^)e" ^( ^ + ^ )7/3^ ' . (33) 



Finally, after rescaling the z' coordinate [z' — > -z'), renaming integration variables, 
and using kinematical identities (^Tj) and Q3"2]), we obtain 



($V)I$M> = J— ^— / dW (0)t (x)</> (0) (x)e-^ , (34) 
V to + 1 7 



where 



This formula was first obtained in |1(J for the semileptonic B — > D^*>ev e decays 



(where = and was already used for the calculation of the corresponding 
Isgur-Wise function |i~0|j-[[l~4] . 



To illustrate the use of the results obtained in this section, we consider the non- 
relativistic quark model as one simple example. 



4 Non-relativistic quark model 

Assuming that we can describe heavy-light mesons using a simple non-relativistic 
potential model, the rest frame LDF wave functions (with angular momentum j and 
its projection Xj), can be written as 

fijl? ( x ) = S - R «i( r ) 1 imL( fi )Xm s (-^,mL;-,m s |j,A i ;L, -) , (36) 

m L ,m a * * 
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where Xm a represent the rest frame spinors normalized to one, xL'Xm s = 5 m >^ mg , 
and a represents all other quantum numbers. 

Explicitly, taking into account Clebsch-Gordan coefficients, the states that we 
need are 



2 2 



x) = R a0 Y oXi , (37) 



,(<*!) < ^ - r / 2, „ / .1 



'11 

2 2 



x) = R a x[^Y llX _ h - J -Y loXh ] , (38) 



^(x) 

2 2 



^iJ^iiX-i + J^ioX|] , (39) 



><<?> 

2 2 



x) = i^J|y aiX _ i -./|y aoX4 ]. (40) 



5 U 5 



Now, using the well known expression 



OC 



47T 



e"^ = £(2Z + l)H)'a(*0y , (41) 

together with the wave functions given above, the overlap expression (0) givesQ 



($f?(^)KiM) = J^rO'o(ar))^, (42) 
22 22 y u; + 1 



*tl(v')m(v)) = 'i/— (ii(ar));r . (13) 

($l7(^)l$?iM) = -i^J^-iUar))^ , (44) 
22 22 y a> + 1 



where 



22 22 V OJ+1 



v*V)>2a = / r 2 dr< a ,(r) J R Qi (r)F(r) . (46) 



1 Using the states analogous to (|37|)-([40|), but for Xj = — |, one can verify that all overlaps 
indeed satisfy ( |T3| ) and ([l4|), so that everything is consistent. 
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At this point, one can readily obtain the values for the IW functions (p3|)-(p6|) 
and their derivatives at the zero recoil point. Let us summarize the results (ignoring 
irrelevant phase factors and suppressing quantum numbers a' and a): 

• 07 — > (07,17) transitions. 

2 22 

fcH = -^-rOo(ar))oo , (47) 
u; -(- l 

fc(l) = (l)oo, (48) 

e'c(i) = -^-^(^- + 4) 2 <- 2 >oo • (49) 

Note that these expressions include transitions from the ground state into 
radially excited states. If the two j = \ states are the same, £c(l) is normalized 
to one and E'g = Eg. 

• 07 — > (Ot, it) transitions. 

2 22 

( E (u) = J_ - Oi(ar)) M , (50) 
— 1 

&(1) = ^- + 4)(r> 10 , (51) 
^(1) = -^(i? 5 + 4)(r)io-^(i? 9 - + 4) 3 <r 3 > 1 o . (52) 



0i — ► (13,23) transitions. 

2 22 



/ 3 2 

C-f(^) = J 2 _ luj+ 1 0'iK))io , (53) 

Cf(1) = ^(^- + 4)(r) 10 , (54) 

et-(i) = -^(^-+4)( r )io-^(^-+4) 3 <^ 3 >io • (55) 
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Oj — > (1 3 ,2 3 ) transitions. 



-3)^3 
2 2 2 



&?M = (j2(ar)) 20 , (56) 

ur — 1 



fo(l) = I ^ f (^- + 4) 2 (r 2 ) 20 , (57) 
= - T ^(^- + ^ 2 (r 2 ) 20 -^^(i? 5 + 4) 4 <r 4 > 2 o • (58) 

Let us briefly compare these results with the formulation of ||, which was also 
followed by [[J. There, apart from irrelevant phase factors, different IW functions 
are identified directly as overlaps of the wave functions of the initial and the final 
state of the LDF in the rest frame of the initial meson. Explicitly (putting a tilde 
over the form factors to avoid confusion), 

ic(u) = Oo(ar)) 00 , (59) 

£ E (u) = V3(j 1 (ar)) l0 , (60) 

If(u) = V3(jx(ar)) l0 , (61) 

£ G (u) = V5(j 2 (ar)) 20 , (62) 

with the definition 



a = E' q Vuj 2 - 1 . (63) 

This formulation has several difficulties. For example, if one uses harmonic oscillator 
wave functions (as was done in ||) in order to estimate overlaps (jj(ar)), then it is 
easy to show that 



= S ' (64) 



where (5b is the variational parameter. Using the values from the ISGW model |T5| 
(which is also used in ||), (3b ~ 0.4 GeV and E'^ w 330 MeV, one gets the slope of 
the IW function for the semileptonic B decays, 

l' c {\) w -0.34 , (65) 
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which is clearly too large [jj"E| . On the other hand, the same wave function used in 
(g^) gives (Eg = E' ? « 330 MeV) 

C c (l) « -0.84 , (66) 

a result which is in much better agreement with the data [fL6| . Furthermore, from 
|62|) one can see that all form factors (except for £c) vanish at the zero recoil 



point, whereas the Bjorken sum rule |T7| requires a nonvanishing P-wave form factor 
in this limit. It is clear that our form factors fl50|) and (|53f) do not suffer from that 
problem. Also, from (|60|) and (|6l|) one can see that £e(u;) = £f(^)) while ([5(]) and 



imply that our form factors satisfy 

= , (67) 

and in particular 

e*(i) = ^fr(i) ■ (68) 

Even though we shall postpone a full account of the radiative rare B decays for later 



18], let us just point out that the authors of pj emphasize a substantially larger 



branching fraction for the decay B — > (1430)7 than the one found in ||19|| . If one 



includes only the factor of y~Y^Y from (|25| ) (which is about ^= for to ph 2), one 
gets a factor of 3 smaller result than the one quoted in ||, bringing this particular 



branching ratio into much better agreement with the result of [19 . 

Finally, we can generalize the quark model approach to any model involving the 
Dirac equation with a spherically symmetric potential. There, the wave function 
has the form 



(69) 



where y^ x . are the usual spherical spinors, k — I (I — j + |) or k — — I — 1 (I — j — |), 
and a again denotes all other quantum numbers. Using this it can be shown that all 
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the expressions p^)-p5|) and (^7|)-([58|) are unchanged, except that the expectation 
value is replaced by 



(F(r))il - (F(r))ff = J r 2 dr[rj;,(r)f^(r) + g*J f (r)g k aj (r)]F(r) . (70) 

5 Conclusions 

We have presented a simple method for relating form factors as defined by the 
covariant trace formalism || |6], to the explicit overlaps of the rest frame wave 
functions describing the initial and the final states of the light degrees of freedom. 
We have obtained explicit formulae for several transitions of interest (from 0~ into 
a few lowest excited states), and have shown how one can apply these expressions in 
the simple quark model, and in models involving the Dirac equation with spherically 
symmetric potentials. We have also pointed out several inconsistencies present in 
recent calculations of radiative rare B decays into higher K- resonances [|[ ||, and 
have shown how these can be fixed even within the same non-relativistic quark model 
that was used in H. A full account of the radiative rare B decays will be presented 



elsewhere 18 . 
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